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Numerical  Analysis  of  Third  Harmonic  Overtone  of 
Thickness-Shear  Vibrations  in  SC-Cut  Quartz  Resonators. 

Final  Report 


Yook-Kong  Yong 

Dept,  of  Civil  and  Environmental  Engineering, 
Rutgers  University,  Piscataway,  New  Jersey. 


I.  Statement  of  the  problem  studied. 

Research  and  development  of  numerical  methods  and  visualization  techniques  for  the 
analysis  high  frequency  SC-cut  quartz  resonators  vibrating  at  the  third  harmonic  overtone 
of  thickness-shear  mode. 


II.  Summary  of  most  important  results. 

1.  Detailed  numerical  modeling  and  visualization  of  high  frequency  modes  of  vibration  in 
quartz  crystal  resonators  such  as  that  of  the  fundamental  thickness  shear  mode  and  its 
harmonic  overtones  were  implemented.  Software  and  hardwares  for  visualization  of 
very  high  frequency  plate  vibrations  were  investigated,  \dsualization  techniques  for  2- 
D,  2.5-D  and  3-D  data  sets  were  explored.  The  development  and  results  were  reported 
in  publications  1  to  5. 

2.  Theoretical  study  and  development  of  two-dimensional  piezoelectric  plate  equations 
for  the  numerical  modeling  of  very  high  frequency  piezoelectric  crystal  resonators  were 
completed.  Dispersion  cmves  for  an  infinite  piezoelectric  plate  with  friee  surfaces  were 
calculated  using  both  the  three-dimensional  piezoelectric  equations,  and  the  two- 
dimensional  piezoelectric  plate  equations,  and  compared.  'Die  dispersion  curves  were 
calculated  up  to  the  third  overtone  of  thickness  shear.  A  good  fit  between  the  dispersion 
relations  obtained  from  the  two-dimensional  piezoelectric  plate  equations  and  the 
three-dimensional  theory  indicates  that  the  two-dimensional  equations  will  yield 
accmate  results.  The  comparison  showed  a  relatively  good  fit  for  the  real  wave 
numbers  which  represented  the  propagating  modes.  The  fit  for  the  imaginary  wave 
numbers  which  represented  the  evanescent  modes  were  not  as  good.  The  third  overtone 
SC-cut  quartz  resonator  will  vibrate  in  the  vicinity  of  a  normalized  frequency  equal  to 
2.8.  Straight  crested  wave  studies  of  SC-cut  quartz  plate  resonators  vibrating  at  the 
third  harmonic  overtone  of  thickness-shear  mode  were  completed.  The  plate  model 
incorporated  all  mechanical  displacement  and  electric  potential  components  up  to  the 
third  order  of  Lee,  Syngellakis  and  Hou’s  piezoelectric  plate  equations.  Hence,  there 
were  12  mechanical  displacement  components,  and  four  electric  potential  components. 
The  mass  and  electrical  effects  of  a  pair  of  electrode  conducting  films  on  the  top  and 


bottom  surfaces  of  the  plate  resonator  were  studied.  The  use  of  Muller’s  method  in  the 
calculation  of  the  dispersion  relations  is  critical  if  there  were  more  than  six  components 
of  displacements  in  the  plate  model.  The  results  were  reported  in  publication  7. 

3.  Finite  element  implementation  of  the  piezoelectric  plate  equations  and  numerical  ^ 
modeling  of  SC-cut  quartz  resonators  vibrating  at  the  third  overtone  of  thickness  shear 
mode  were  performed.  The  results  were  reported  in  publications  8  and  13. 

4.  From  our  working  experience  on  the  finite  element  modeling  of  high  frequency  SC-cut 
quartz  resonators,  we  learned  some  important  strategies  and  problems  in  such 
numerical  models: 

(a)  For  accurate  representation  of  eigenpairs,  the  number  of 
elements  per  half  wave  of  eigenmode  must  be  at  least  three  or  more. 

(b)  The  Lanczos  algorithm  for  the  eigenvalue  problem  was  found  to 
be  superior  in  terms  of  accuracy  and  computational  speed.  The 
algorithm  allowed  for  the  calculation  for  a  few  eigenpairs  in  a 
narrow  bandwidth  at  any  frequency. 

(b)  For  the  piezoelectric  eigenvalue  matrix  problem,  better 
algorithms  were  needed  for  managing  storage  memory  and 
increasing  computational  speed.  This  was  the  motivation  for  the 
next  phase  of  research. 

5.  Research  on  new  algorithms  for  the  solution  of  the  piezoelectric  eigenvalue  matrix 
problem  were  performed.  The  customary  method  was  to  statically  condense  the 
dielectric  matrix  since  the  problem  is  electrostatic.  The  static  condensation  however 
destroyed  the  sparse  matrix  stracture  of  the  stiffness  matrix,  and  caused  unacceptable 
increase  in  memory  storage  and  computational  time.  Two  schemes  were  proposed: 

(a)  A  perturbation  scheme  in  which  the  mechanical  eigenvalue 
problem  was  solved  independently  and  the  piezoelectric  effects 
were  subsequently  added  as  a  perturbation.  This  scheme  was 
however  applicable  only  to  materials  with  very  weak  piezoelectric 
coupling  such  as  quartz.  It  would  not  be  applicable  to  materials 
such  as  lithium  niobate,  zinc  oxide,  and  aluminum  nitride.  The 
results  were  reported  in  publication  10.  A  new  method  for  storing 
the  mass  matrix  which  reduced  the  memory  storage  by  about  94% 
was  also  presented 

(b)  A  scheme  which  avoided  static  condensation  of  the  dielectric 
matrix  was  implemented.  This  scheme  was  superior  to  the 
perturbation  scheme  because  it  was  applicable  to  all  piezoelectric 
materials,  was  simpler  conceptually,  and  was  more  efficient  in 
terms  of  memory  and  computational  speed.  The  old  customary 
scheme  of  storing  the  electrical  degrees  of  freedom  separate  from 
the  mechanical  degrees  of  freedom  was  discarded.  A  new  storage 
scheme  which  mixed  the  mechanical  and  electrical  degrees  of 
freedom  in  the  global  matrix  so  as  to  reduce  its  half-bandwidth  was 
implemented.  This  resulted  in  savings  of  more  than  two  orders  of 
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magnitude  in  memory  and  computational  speed.  Part  of  the  results 
were  reported  in  publication  14.  Further  results  will  be  reported  in 
the  Proceedings  of  the  IEEE  1995  Ultrasonics  Symposium. 

6.  Implementation  of  the  finite  element  method  for  piezoelectric  laminated  plate  witl^^ 
electrodes  was  performed.  The  studies  were  reported  in  publications  6, 9, 1 1, 12  and  15. 

7.  The  accuracy  of  plate  equations  for  straight  crested  waves  at  frequencies  in  the  vicinity 
of  the  fundamental  thickness  shear  mode  in  an  SC-cut  quartz  strip,  was  examined.  The 
first  order  plate  theories  with  correction  factors  were  currently  assumed  to  predict 
accurately  the  plate  modes  of  vibration  at  frequencies  up  to  twenty  percent  higher  than 
the  fundamental  thickness  shear  frequency.  These  plate  equations  were  derived  using  a 
series  expansion  of  the  thickness  modes  with  various  functions  of  the  thickness 
coordinate:  (1)  Power  series  (Mindlin),  (2)  trigonometric  series  (Lee,  et.  al.),  (3) 
Legendre  series  (Mindlin,  et.  al.)  and  (4)  normal  modes  (Peach).  Our  results  showed 
that  the  first  order  plate  equations  did  not  yield  accurate  frequency  spectra  of  the  modes 
in  the  vicinity  of  the  fundamental  thickness  shear  mode,  although  the  thickness  shear 
mode  itself  was  predicted  accurately.  Four  finite  element  programs  were  used  in  the 
study.  The  programs  individually  calculated  the  straight  crested  wave  frequency  spectra 
as  a  function  of  the  length  to  thickness  ratio  for  an  SC-cut  quartz  strip  using  (1)  the  first 
order  Mindlin  plate  equations,  (2)  the  third  order  Mindlin  plate  equations,  (3)  the  third 
order  Lee,  et  al.  plate  equations,  and  (4)  the  three  dimensional  elastic  equations.  The 
numerical  solutions  for  the  three  dimensional  elastic  equations  were  checked  against 
the  simple  analytical  solutions  for  an  elastic  membrane.  The  results  showed  that  all  the 
plate  theories  predicted  the  fundamental  thickness  mode  accurately,  but  the  frequency 
spectra  of  modes  in  the  vicinity  are  not  predicted  well  by  the  first  order  Mindlin  plate 
equations  and  the  third  order  Lee  et.  al.  plate  equations.  The  third  order  Mindlin  plate 
equations,  on  the  other  hand,  predicted  very  well  the  entire  frequency  spectrum  in  the 
vicinity.  Hence,  if  an  accurate  prediction  of  the  unwanted  modes  near  the  fundamental 
thickness  shear  mode  was  needed,  the  third  order  Mindlin  plate  equations  were 
recommended.  The  third  order  Mindlin  plate  equations  did  not  require  shear 
corrections  factors.  The  mode  shapes  of  the  third  order  displacement  components 
showed  a  surprisingly  large  contribution  to  the  resonant  modes  in  the  neighborhood  of 
the  fundamental  thickness  shear  frequency.  The  results  will  be  reported  in  the 
Proceedings  of  the  IEEE  1995  Annual  Symposium  on  Frequency  Control. 


III.  Interactions  with  the  U.S.  Army  Research  Laboratory, 
Electronics  &  Power  Sources  Directorate,  Fort  Monmouth. 

Interactions  are  maintained  with  Dr.  Arthur  Ballato  and  Dr.  John  R.  Vig  of  the  U.S.  Army 
Research  Laboratory,  Electronics  &  Power  Sources  Directorate,  Fort  Monmouth.  There 
are  current  works  on  acceleration  sensitivity  of  SC-cut  quartz  crystal  resonators  at  the  said 
laboratory.  Numerical  modeling  and  visualization  of  very  high  frequency  quartz  crystal 
resonators  can  be  applied  to  acceleration  sensitivity  studies. 
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Visualization  Techniques  for  Bulk  Wave  Resonators 


D.  Silver  J.  T.  Stewart 

Y-K.  Yong 
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Rutgers  University 
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Abstract 

In  this  paper,  we  present  an  overview  of  some  fun¬ 
damental  algorithms  for  visualizing  multi- dimensional 
data  sets.  These  techniques  are  applied  to  ongoing  re¬ 
search  problems  in  bulk  wave  piezoelectric  resonators. 

Introduction 

The  goal  of  large  scale  simulations  and  experiments 
in  computational  science  is  a  quantitative  and  math¬ 
ematical  understanding  of  the  model  being  investi¬ 
gated.  A  key  component  in  this  process  is  visual¬ 
ization:  the  use  of  computer  graphics  to  depict  and 
project  the  resulting  data.  As  experiments  become 
more  sophisticated,  the  data  generated  becomes  over¬ 
whelming  and  new  techniques  are  needed  to  help  un¬ 
derstand  and  aid  in  its  interpretation. 

The  sequence  of  steps  involved  in  this  process  can 
be  summarized  as  follows: 

•  Data  generation  (simulation,  etc.) 

•  Filtering,  i.e.  deriving  other  variables,  interpolci- 
tion,  etc. 

♦  Creating  geometric  primitives:  converting  the 
data  to  a  form  suitable  for  rendering, 

#  Rendering:  Displaying  the  geometry  (with  shad¬ 
ing,  transparency...) 

This  is  process  is  not  a  static  one,  i.e.  the  process 
may  involve  many  iterations  before  a  suitable  visual¬ 
ization  results. 

Scientific  visualization  is  an  active  and  growing 
field.  Almost  all  domains  of  science  are  incorporat¬ 
ing  computer  graphics  algorithms  into  standard  scien¬ 
tific  experiments.  These  fields  include  computational 
fluid  dynamics,  medical  imaging,  molecular  biology, 
chemistry,  meteorology,  etc.  Although  each  scientific 
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domain  has  specific  requirements,  most  of  the  proce¬ 
dures  used  are  generic  and  can  be  adapted  to  other 
fields.  Now  that  finite  element  solutions  (2D  and  3D) 
are  more  common  for  bulk  wave  resonators,  visualiza¬ 
tion  techniques  are  becoming  increasingly  important 
to  help  in  the  analysis  of  the  resulting  data  sets. 

In  this  paper,  we  present  an  overview  and  imple¬ 
mentation  details  of  some  fundamental  algorithms  for 
visualizing  multi-dimensional  data  sets.  These  tech¬ 
niques  are  applied  to  ongoing  research  problems  in 
bulk  wave  piezoelectric  resonators.  These  methods 
may  be  useful  for 

♦  Identifying  modes:  for  example,  thickness  shear, 
overtones  (Figure  2),  and  energy  trapping  phe¬ 
nomena  (Figure  3); 

♦  Detecting  symmetry:  whether  a  mode  is  symmet¬ 
ric  or  anti-symmetric  (Figure  13); 

♦  Determining  mounting  strategics:  locate  zero- 
values  (Figure  7). 

In  the  first  section,  we  present  a  set  of  two  di¬ 
mensional  routines.  For  modeling  a  vibrating  res¬ 
onator  using  3D  finite  elements,  a  combination  of  slic¬ 
ing,  isosurface  contours,  and  volume  visualization  are 
used.  These  techniques  arc  also  applicable  to  layer- 
wise  plate  theory.  We  hope  to  show  that  using  visual¬ 
ization  tools  can  greatly  enhance  the  process  of  scien¬ 
tific  discovery  and  progress  in  this  area  of  research. 

(Note:  The  pictures  presented  in  this  paper  were 
originally  in  color  -  some  of  the  effects,  especially 
transparency,  do  not  reproduce  well  in  black  and 
white.) 

Two  Dimensional  Techniques 

This  section  discusses  methods  for  visualizing  2D 
scalar  fields  (i.e.  an  x,y  grid  with  a  value  at  each 
node)  resulting  from  finite  clement  solutions  for  the 
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Mass-Frequency  Influence  Surface,  Mode  Shapes, 
and  Frequency  Spectrum  of  a  Rectangular 
^T-Cut  Quartz  Plate 

Yook-Kong  Yong,  Member,  IEEE  and  Janies  T.  Stewart 


Abstract— The  mass-frequency  influence  surface  and  frequency  spec¬ 
trum  of  a  rectangular  /IT-cut  quartz  plate  is  studied.  The  mass-fre- 
queiicy  influence  surface  is  defined  as  a  surface  giving  the  frequency 
change  due  to  a  small  localized  mass  applied  on  the  plate  surface.  Fi¬ 
nite-element  solutions  of  Mindlin’s  two-dimensional  (2-1))  plate  equa¬ 
tions  for  thickness-shear,  thickness-twist,  and  flexural  vibrations  are 
given.  Spectrum  splicing,  and  an  efficient  eigenvalue  solver  using  the 
Lanezos  algorithm  were  incorporated  into  the  finite-element  program, 

A  convergence  study  of  the  fundamental  thickness-shear  mode  and  its 
first  symmetric,  anhannonic  overtone  was  performed  for  finite-ele¬ 
ment  meshes  of  increasing  fineness.  As  a  general  rule,  more  than  two 
elements  must  span  any  half-wave  in  the  plate  or  spurious  mode  shapes 
will  be  obtained.  Two-dimensional  (2-1))  mode  shapes  and  frequency 
spectrum  of  a  rectangular  AT-cut  plate  in  the  region  of  the  fundamen¬ 
tal  thickness-shear  frequency  are  presented.  The  mass-frequency  influ¬ 
ence  surface  for  a  5-MHz  rectangular,  AT-cut  plate  with  patch 
electrodes  is  obtained  by  calculating  the  frequency  change  due  to  a 
small  mass  layer  moving  over  the  plate  surface.  The  frequency  change 
is  proportional  to  the  ratio  of  mass  loading  to  mass  of  plate  per  unit 
area,  and  is  confined  mostly  within  the  electrode  area,  where  the  mag¬ 
nitude  is  of  the  order  10^  Hz/g. 

I.  Introduction 

HE  PRIMARY  EFFECT  of  mass  loading  on  a  thickness- 
shear  resonator  is  well  known,  that  is,  the  resonator  fre¬ 
quency  is  lowered  in  proportion  to  the  magnitude  of  mass  load¬ 
ing.  It  was  shown  theoretically  [1],  [21  that  very  small, 
fluctuating  mass  loading,  namely  that  of  a  fluctuating  mono- 
layer  of  molecules,  can  have  a  measurable  effect  on  the  phase 
fluctuations  of  ultra-high-frequency  resonators.  In  this  paper, 
the  frequency  effect  of  small  mass  loadings  on  a  5-MHz,  rect¬ 
angular  /IT-cut  quartz  plate  with  rectangular  electrode  films  is 
calculated  using  finite  clement  solutions  of  Mindlin’s  plate 
equations  with  mass  loadings  [3].  The  plate  resonator  is  vibrat¬ 
ing  in  the  fundamental  thickness-shear  mode.  The  frequency 
spectrum  and  mode  shapes  of  resonant  modes  in  the  region  of 
the  fundamental  thickness-shear  mode  are  also  presented. 

11.  Mindlin’s  Plate  Equations  for  Thickness- 
Shear,  Thickness-Twist,  and  Flexural 
Vibrations  in  an  AT-cut,  Electrode 
Plated  Quartz  Plate 

The  calculation  of  frequency  change  due  to  a  small,  discrete 
mass  loading  on  a  plate  surface  lends  itself  naturally  to  the  finite 
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element  method.  Mindlin’s  equations  for  an  AT-cut,  electrode 
plated  quartz  plate  [3],  14]  are  employed.  Lee,  Zee,  and  Brcbbia 
[5]  employed  these  equations  in  their  study  of  the  influence  of 
the  size  and  mass  of  rectangular  patch  electrodes  on  the  vibra¬ 
tions  and  energy  trapping  in  a  rectangular  AT-cut  quartz  plate. 
Fig.  1  shows  a  rectangular  crystal  plate  of  density  p,  thickness 
2b,  length  21,  and  width  2w  that  is  partially  plated  on  the  top 
and  bottom  surfaces  by  electrode  platings  of  density  p',  thick¬ 
ness  2b\  length  2a,  and  width  2c  over  the  central  portion  of  its 
faces. 

A.  Stress  Equations  of  Motion 

The  stress  equations  of  motion  accommodating  thickness- 
shear,  thickness-twist,  and  flexural  vibrations  arc 

Qu  +  Qi.3  =  2hpil  +  R)ii„  (1) 

M,.,  +  M5.3  -  2i  =  5  (2) 

W5.1  +  A^3.3  -  23  =  5  +  2R)h,  (3) 

where  the  terms  Q,  (/  =  1,  3)  are  the  transverse  shear  forces 
per  unit  length,  and  M,-,  for  f  =  1,  or  3,  are  the  bending  mo¬ 
ments  per  unit  length,  and,  M5  is  the  twisting  moment  per  unit 
length.  The  term  R  =  2p'b'  f(pb)  is  ratio  of  mass  of  electrodes 
to  mass  of  plate  per  unit  area.  Displacement  U2  is  the  transverse 
displacement  in  the  X2  direction,  and  i/^i,  and  1/^3  are  shear  ro¬ 
tations  about  the  jcj,  and  axes,  respectively. 

B.  Stress-Strain  Relations 

For  an  AT-cut  quartz  plate,  the  clastic  constants  exhibits 
monoclinic  symmetry  when  the  x^  axis  coincides  with  one  of 
the  digonal  axes  of  crystal.  The  constitutive  relations  arc 


(2i  —  2bk]C(y(,yi 

H) 

03  ~  2M3C4473 

(5) 

M,  =  lb\CuX^  +  C,3X3) 

(6) 

A^3  ”  lb^{^\3X\  +  Q3X3) 

(7) 

M5  “  \b^Cf,sX5  , 

(8) 

where  7/  (/  ==  1,  3)  are  the  shear  strains,  and  Xn  X31  X5 
are  thc.t,  bending  curvature,  .^3  bending  cimvaturc^and  twisting 
curvature,  respectively.  The  cocflicicnls  k]  and  k]  in  (4)  and 
(5)  are  shear  correction  factors  tliat  include  mass  loading  cfl'ccts 
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Abstract — Finite-element  solutions  for  the  fundamental  thick¬ 
ness  shear  mode  and  the  second-anharmonic  overtone  of  a  cir¬ 
cular,  1.87  MHz  AT-cut  quartz  plate  with  no  electrodes  are 
presented  and  compared  with  previously  obtained  results  for  a 
rectangular  plate  of  similar  properties.  The  edge  flexural  mode 
in  circular  plates,  a  vibration  mode  not  seen  in  the  rectangular 
plate,  is  also  presented.  A  5-MHz  circular  and  electroded  AT-cut 
quartz  plates  is  studied.  A  portion  of  the  frequency  spectrum  is 
constructed  in  the  neighborhood  of  the  fundamental  thickness- 
shear  mode.  A  convergence  study  is  also  presented  for  the 
electroded  5-MHz  plate.  A  new  two-dimensional  (2-D)  technique 
for  visualizing  the  vibration  mode  solutions  is  presented.  This 
method  departs  substantially  from  the  three-dimensional  (3-D) 
“wire-frame”  plots  presented  in  the  previous  analysis.  The  2-D 
images  can  be  manipulated  to  produce  nodal  line  diagrams  and 
can  be  color  coded  to  illustrate  mode  shapes  and  energy  trapping 
phenomenon.  A  contour  plot  of  the  mass-frequency  influence 
surface  for  the  plated  5-MHz  resonator  is  presented.  The  mass- 
frequency  influence  surface  is  defined  as  a  surface  giving  the 
frequency  change  due  to  a  small  localized  mass  applied  to  the 
resonator  surface. 


I.  Introduction 

The  use  of  quartz  crystal  plates  in  ultrasonics  is  wcll^ 
studied  experimentally.  Analytical  closed  form  solutions 
of  a  freely  vibrating  anisotropic  elastic  plate  can  in  general 
only  be  obtained  for  either  a  one-dimensional  (1-D)  problem, 
or  in  an  infinite  domain.  For  a  two-dimensional  (2-D)  finite 
domain,  solutions  may  be  obtained  either  approximately,  or  for 
certain  specific  types  of  boundary  conditions  and  geometries. 
To  perform  a  practical  analysis  of  the  AT-cut  quartz  resonator, 
it  is  necessary  to  solve  the  free  vibration  problem  for  a  finite 
plate.  In  this  paper,  such  a  problem  is  solved  using  the  finite- 
clement  method. 

The  high  frequency  vibrations  as  well  as  the  anisotropicity 
of  quartz  demands  that  a  general  thick  plate  theory  be  used. 
In  this  study,  Mindlin’s  2-D  plate  equations  are  truncated 
to  a  first-order  approximation.  The  first-order  approximation 


X. 


Fig.  1.  Oricntalion  and  dimensions  of  a  circular  electroded  plate. 


Fig.  2.  Finite-element  mesh  of  quarter  circular  plate  with  essential  boundary 
conditions  at  the  plate  edges  along  the  coordinate  axes. 
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is  chosen  because  il  is  the  lowest-order  formulation  which 
includes  the  major  modes  of  interest,  namely,  the  fundamental 
thickness-shear  modes  [2],  Real  /f'/'-cul  resonators  generally 
contain  a  thin  electrode  plating  on  a  portion  of  each  face, 
across  which  a  time  varying  voltage  is  applied  to  drive  the 
vibration.  The  effect  of  this  electrode  plating  was  studied 
analytically  by  Mindlin  [3]  for  an  infinite  plate.  Mindlin’s 
equations  for  the  electroded  plate  are  incorporated  into  the 
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Abstract 

Finite  element  matrix  equations  employing 
high  frequency,  piezoelectric  plate  equations  are 
derived.  The  equations  may  be  used  for  modeling 
third  harmonic  overtone  of  thickness-shear 
vibrations.  A  perturbation  technique  is  developed 
to  account  for  piezoelectric  stiffening  in  the 
mechanical  stiffness  matrix.  Results  from  the 
perturbation  method  compares  well  with  the 
direct  solution  of  the  piezoelectric  finite  element 
equations.  The  technique  will  result  in  significant 
.savings  in  computer  memory  and  computational 
lime.  Numerical  results  for  straight  crested 
waves  in  a  third  overtone  SC-cut  quartz  strip 
with  and  without  electrodes  are  presented. 

Steady  slate rcspon.se  to  an  electrical  excitation  is 
calculated. 

1.  Introduction 

We  discuss  in  this  paper  the  finite  element 
fonnulation  and  results  on  straight  crested  waves 
in  the  vicinity  of  the  third  harmonic  overtone  of 
thickness  shear,  piezoelectric  vibrations  in  a  SC- 
cut  quartz  strip.  The  straight  crested  waves  are 
spatially  one-dimensional,  but  since  the  study  is 
])art  of  it  wider  study  of  third  overtone  SC-cut 
quartz  resonators,  two-dimensional  finite 
clement  equations  for  piezoelectric  plate 
vibrations  are  derived  and  presented.  A  one- 
dimensional  finite  element  program  is  developed 
which  could  be  modified  subsequently  to  handle 
two-dimensional  problems.  The  mechanical 


stiffness  of  the  finite  element  matrix  equations 
for  the  SC-cut  quartz  strip  is  a  few  orders  of 
magnitude  larger  than  the  piezoelectric  stiffness, 
and  the  piezoelectric  coupling  is  weak.  Hence, 
the  piezoelectric  effects  could  be  modelled  using 
a  perturbation  technique.  A  new  perturbation 
method  for  piezoelectric  vibrations  is  derived 
and  proposed. 

II.  Plate  Equations  for  High  Frequency 
Piezoelectric  Vibrations. 

The  piezoelectric  plate  equations  for  high 
frequency  vibrations  employed  in  the  cunent 
finite  element  formulation  were  derived  by  Lee, 
Syngellakis  and  Hou  [1].  The  components  of 
mechanical  displacement  and  electric  potential 
are  expanded  in  an  infinite  series  with  their 
thickness-dependence  expressed  by 
trigonometric  functions: 

M.  UpXj.Xj, /)  =  (XpXj,  0COS^(l  -^) 

(1 

0  =  Ul.Xj,  OCOSy  (1  -^) 

n  =  0 

X2 

where  ^  =  -r  is  a  nondimensional  thickness 
0 

coordinate  of  the  plate.  The  terms  uj"^  and 
are  functions  of  spatial  coordinates  X]  and  X3 
only  and  time  t.  We  present  without  formal 
derivations  the  following  field  equations  from 
reference  [1].  The  interested  reader  may  obtain 
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Numerical  Algorithms  and  Results  for 
SC-Cut  Quartz  Plates  Vibrating  at  the  Third 
Harmonic  Overtone  of  Thickness  Shear 
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Abstract — Finite  element  matrix  equations,  derived  from  two- 
dimensional  piezoelectric  high  frequency  plate  theory  are  solved 
to  study  the  vibrational  behavior  of  the  third  overtone  of  thick¬ 
ness  shear  in  square  and  circular  SC-cut  quartz  resonators.  The 
mass-loading  and  electric  effects  of  electrodes  are  included.  A 
perturbation  method  which  reduces  the  memory  requirements 
and  computational  time  significantly  [2]  is  employed  to  calculate 
the  piezoelectric  resonant  frequencies.  A  new  storage  scheme 
is  introduced  which  reduces  memory  requirements  for  mass 
matrix  by  about  90%  over  that  of  the  envelope  storage  scheme. 
Substructure  techniques  are  used  in  eigenvalue  calculation  to 
save  storage.  Resonant  frequency  and  the  mode  shapes  of  the 
harmonic  third  overtone  thickness  shear  vibrations  for  square 
and  circular  plates  are  calculated.  A  predominant  third  overtone 
thickness  shear  displacement,  coupled  with  the  third  overtone  of 
thickness  stretch  and  thickness  twist,  is  observed.  Weak  coupling 
between  the  third  order  thickness  shear  displacement  and  the 
zeroth-,  first-,  and  second-order  displacements  is  noted.  The 
magnitudes  of  the  lower  order  displacements  are  found  to  be 
about  two  orders  smaller  than  that  of  the  third  overtone  thickness 
shear  displacement. 

L  Introduction 

This  paper  describes  computational  strategies  in  the  finite 
element  analysis  of  SC-cut  quartz  resonators,  and  presents 
the  numerical  results  of  the  third  overtone  thickness  shear 
modes  in  SC-cut  quartz  plates.  Theoretical  developments  and 
applications  of  the  high  frequency  piezoelectric  plate  equations 
were  presented  in  [1]  and  [2],  and  would  not  be  repeated  here. 
The  interested  reader  should  refer  to  the  two  references  for 
further  details.  One  obstacle  in  the  finite  element  modeling  of 
high  frequency  plate  vibrations  is  the  limitation  of  computer 
memory.  A  finite  element  model  using  two-dimensional  (2-D) 
piezoelectric  plate  equations  provides  an  excellent  approxima¬ 
tion  of  the  vibrational  behavior  of  plate  resonators,  and  at  the 
same  time  requires  less  memory  than  the  three-dimensional 
finite  element  model.  The  characteristics  of  a  good  finite 
element  program  rest  on  its  memory  saving  strategies  as  well 
as  its  computational  efficiency. 

The  finite  element  matrix  equations  are  based  on  the  piezo¬ 
electric  plate  theory  derived  by  Lee,  Syngellakis,  and  Hou  [1], 
The  plate  theory  is  chosen  for  its  accuracy  in  predicting  cutoff 
frequencies  of  harmonic  overtones  of  thickness  vibrations. 
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Twelve  mechanical  displacement  and  four  electric  potential 
components,  and  {j  =  1,2,3,  and  n  =  0, 1,2,3), 
are  incorporated  into  the  numerical  model  for  the  analysis  of 
vibrations  with  frequencies  in  the  vicinity  of  the  harmonic 
third  overtone  thickness  modes.  The  mass-loading  and  electric 
effects  of  platings  are  also  included.  In  piezoelectric  vibrations, 
the  terms  in  the  stiffness  matrix  consist  of  two  parts,  namely, 
the  mechanical  stiffness  and  piezoelectric  coupling  stiffness. 
The  latter,  even  though  it  is  a  few  orders  of  magnitude  less  than 
the  mechanical  terms  for  SC-cut  quartz,  destroys  the  sparse 
structure  of  the  stiffness  matrix.  A  perturbation  method  [2]  is 
introduced  to  separately  calculate  the  piezoelectric  frequencies 
from  the  purely  mechanical  frequencies  for  quartz  resonators. 

In  the  calculations  of  purely  mechanical  frequencies,  a  new 
storage  scheme  for  the  mass  matrix  is  developed  which  takes 
advantage  of  its  special  structure  and  reduces  the  memory 
requirements  by  about  90%  over  that  of  the  envelope  storage 
scheme  during  the  Cholesky  factorization,  and  a  substructuring 
scheme  is  employed  to  further  diminish  the  storage  demands. 

Resonant  frequencies  of  the  harmonic  third  overtone  thick¬ 
ness  shear  and  the  associated  mode  shapes  are  calculated 
for  rectangular  and  circular  SC-cul  quartz  plate  resonators. 
A  clean  third  overtone  thickness  shear  displacement,  coupled 
with  the  third  overtone  of  thickness  stretch  and  thickness  twist, 
is  observed.  The  magnitudes  of  zeroth,  first  and  second  order 
displacement  components  are  one  to  two  orders  smaller  than 
those  of  the  third-order  displacement  components. 

II.  Numerical  Strategies  for  SC-Cut  Plates  Vibrating 
AT  THE  Third  Overtone  of  Thickness  Shear 

The  finite  element  matrix  equations  for  piezoelectric  vibra¬ 
tions  have  the  following  form  [2]: 

Kod-\-P'^rp  +  Md  =  Fs  +  Fc  (1) 

Pd-Rrlj  =  Qs  +  Qc  (2) 

where  matrices  Ko,P,R,  and  M  are,  respectively,  the  ele¬ 
ment  mechanical  stiffness,  piezoelectric  coupling,  dielectric 
stiffness,  and  mass.  The  vectors  Fs,Fc,Qs'  Qc 
element  face  traction,  edge  traction,  face  charge,  and  edge 
charge,  respectively,  [2].  The  general  piezoelectric  eigenvalue 
equations  (3)  and  (4)  may  be  derived  from  (1)  and  (2)  by 
setting  the  right-hand  side  vectors  to  zero  and  assuming  the 
harmonic  motion,  namely,  d  =  de"^*  and  rp  = 

Ko'd  +  P'^ip^uj^Md  (3) 
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Exact  Analysis  of  the  Propagation  of  Acoustic  Waves 
in  Multilayered  Anisotropic  Piezoelectric  Plates 

James  T.  Stewart  and  Yook-Kong  Yong,  Member,  IEEE 


Abstract — Exact  analysis  of  the  propagation  of  acoustic  waves 
in  multilayered  piezoelectric  plates  is  performed  using  the  trans¬ 
fer  matrix  method.  A  general  technique  for  analyzing  layered 
piezoelectric  resonators  under  thickness  and  lateral  field  exci¬ 
tation  Is  presented  and  is  applied  to  the  study  of  zinc  oxide 
on  silicon  thin  film  resonators.  Both  one  and  two-dimensional 
analysis  with  general  material  anisotropy  is  performed,  and  a 
simplified  method  for  incorporating  thin  conducting  electrodes 
on  the  plate’s  free  surfaces  is  presented.  The  general  methodol¬ 
ogy  described  is  summarized  into  efficient  algorithms  to  aid  in 
the  implementation  of  the  procedures  and  some  computational 
aspects  are  discussed.  Results  are  presented  for  cutoff  behavior 
as  well  as  general  dispersion  characteristics  for  two  and  three 
layered  plates. 


L  Introduction 

Exact  analysis  of  the  propagation  of  acoustic  waves  in 
multilayered  piezoelectric  plates  is  performed  using  a 
general  and  powerful  transfer  matrix  approach.  The  transfer 
matrix  method  is  a  very  simple  technique  for  analyzing  wave 
phenomena  in  layered  media.  This  technique,  as  applied  to 
piezoelectric  crystals,  has  its  roots  in  the  one  dimensional 
transmission  line  equivalent  circuit  approach  originally  used 
by  electrical  engineers  to  analyze  stacked  crystal  filters.  The 
most  general  analysis  of  this  type  was  given  by  Ballato, 
Bertoni,  and  Tamir  [2].  The  transfer  matrix  method  is  simply 
the  physical  equivalent  of  the  transmission  line  analysis, 
relying  on  the  principles  of  mechanics,  rather  than  on  electric 
circuit  analysis.  Adler  [3]  has  proposed  a  general  method  for 
constructing  transfer  matrices  for  both  bulk  and  surface  wave 
problems  which  is  based  on  linear  systems  analysis.  Nowotny, 
Benes,  and  Schmid  [4]  developed  a  general  one  dimensional 
transfer  matrix  method  for  calculating  the  admittance  matrix 
for  stacked  piezoelectric  layers  with  arbitrary  electrode  place¬ 
ment,  and  Nayfeh  [5]  used  a  transfer  matrix  approach  to  solve 
two  dimensional  problems  of  wave  propagation  at  arbitrary 
angles  in  multilayered  nonpiezoelectric  plates  with  monoclinic 
(or  higher)  symmetry.  The  present  study  represents  a  general 
one  and  two  dimensional  analysis  of  layered  piezoelectric 
plates  with  general  anisotropy.  The  one  dimensional  analysis 
presented  is  similar  to  the  discussion  of  [4]  with  extensions  to 
include  lateral  field  excitation  as  well  as  a  simplified  technique 
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for  handling  the  effects  of  thin  conducting  electrodes  on  the 
plate’s  free  surfaces.  The  two  dimensional  analysis  presented 
represents,  in  principle,  an  extension  of  [5]  to  the  more  general 
case  of  piezoelectric  layers.  The  two  dimensional  analysis 
is  carried  out  for  the  cases  of  both  open  and  short  circuit 
conditions  on  the  plate’s  free  surfaces. 

By  definition,  a  transfer  matrix  is  a  linear  transformation  that 
maps  a  specified  vector  of  field  quantities  from  one  point  to 
another  in  a  given  material.  When  applied  to  a  layered  plate, 
this  method  greatly  simplifies  the  analysis  by  guaranteeing 
continuity  of  field  quantities  across  layer  interfaces  while  re¬ 
ducing  the  problem  to  one  involving  these  quantities  evaluated 
at  the  plate’s  free  surfaces  only.  More  specifically,  a  transfer 
matrix  for  each  layer  is  developed  which  maps  displacements, 
stress  tractions,  electric  potential,  and  electric  charge  from 
the  layers  lower  surface  to  it’s  upper  surface.  Satisfaction 
of  inter-layer  continuity  of  these  field  quantities  leads  to  a 
simple  matrix  multiplication  rule,  from  which  follows  a  global 
plate  transfer  matrix.  The  plate  transfer  matrix  equation  that 
results  is  represented  in  terms  of  the  field  quantities  evaluated 
on  the  plate’s  free  surfaces  only.  From  this,  many  problems 
involving  thickness  field  excitation  of  layered  piezoelectric 
resonators  may  be  handled  quite  readily.  It  is  shown  that, 
with  some  modifications,  this  method  can  handle  problems 
involving  lateral  field  excitation  as  well.  A  simple  technique 
for  incorporating  the  mechanical  effects  of  thin  conducting 
electrodes  on  the  plates  top  and  bottom  surfaces  is  also 
presented. 

The  present  study  involves  the  application  of  the  transfer 
matrix  method  to  the  exact  analysis  of  simple  thickness  modes 
and  general  dispersion  behavior  in  zinc  oxide  on  Silicon  thin 
film  resonators.  Fundamental  resonant  and  antiresonant  cutoff 
frequencies  through  both  thickness  and  lateral  field  excitation 
are  calculated  for  two  and  three  layered  plates.  Along  with 
these  frequencies,  thickness  mode  shapes  for  displacements, 
stress  tractions,  electric  potential,  and  electric  displacement 
are  shown.  Exact  dispersion  relations  for  propagating  straight 
crested  waves  are  also  presented  for  these  plates.  Both  real 
and  imaginary  branches  of  these  dispersion  relations  are  shown 
and  these  calculations  are  carried  out  for  both  open  circuit  and 
short  circuit  conditions.  Two  different  crystal  tensor  rotations 
of  the  ZnO  layers  are  considered,  giving  rise  to  different  cases 
of  electro-mechanical  coupling. 

II.  General  Problem 

Consider  a  plate  of  thickness  2b  constructed  of  N 
anisotropic  piezoelectric  layers  stacked  normal  to  the  X2 
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Abstract — When  the  piezoelectric  stifTening  matrix  is  added 
to  the  mechanical  stiffness  matrix  of  a  finite  element  model,  its 
sparse  matrix  structure  is  destroyed.  A  direct  consequence  of  this 
loss  in  sparseness  is  the  significant  rise  in  memory  and  compu¬ 
tational  time  requirements  for  the  model.  For  weakly  coupled 
piezoelectric  materials,  the  matrix  sparseness  can  be  retained  by 
a  perturbation  method  which  separates  the  mechanical  eigenvalue 
solution  from  its  piezoelectric  effects.  A  perturbation  and  finite 
element  scheme  for  weakly  coupled  piezoelectric  vibrations  in 
•quartz  plate  resonators  has  been  developed.  Finite  element  matrix 
equations  were  derived  specifically  for  third  overtone  thickness 
shear,  SC-cut  quartz  plate  resonators  with  electrode  platings. 
High  frequency,  piezoelectric  plate  equations,  previously  derived 
by  Lee,  Syngellakis,  and  Hou  [1],  were  employed  in  the  for¬ 
mulation  of  the  finite  element  matrix  equation.  The  equations 
may  be  used  for  modeling  third  harmonic  overtone  of  tUckness 
shear  vibrations.  Results  from  the  perturbation  method  for  SC- 
cut  quartz  plates  compared  well  with  the  direct  solution  of  the 
piezoelectric  finite  element  equations.  This  method  will  result 
in  significant  savings  in  computer  memory  and  computational 
time.  Resonance  and  anti-resonance  frequencies  of  a  certain 
mode  could  be  calculated  easily  by  using  the  same  eigen-pair 
from  the  purely  mechanical  stiffness  matrix.  Numerical  results 
for  straight  crested  waves  in  a  third  overtone  SC-cut  quartz 
strip  with  and  without  electrodes  are  presented.  The  steady-state 
response  to  an  electrical  excitation  is  calculated.  The  accuracy 
of  the  perturbation  method  is  dependent  on  the  magnitude  of 
piezoelectric  coupling  constant  in  the  mode  of  vibration.  For 
example,  an  investigation  on  the  first  few  modes  of  vibration 
in  a  cube  made  of  lithium  niobate,  which  has  high  piezoelectric 
coupling  consunts,  yielded  divergent  results  for  the  perturbation 
method.  Hence,  the  method  should  not  be  used  for  materials  with 
high  piezoelectric  coupling  constants. 


1.  iNTRODUenON 

WE  discuss  in  this  paper  a  perturbation  method  for 
weakly  coupled  piezoelectric  vibrations,  the  finite  ele¬ 
ment  formulation,  and  results  on  straight  crested  waves  in  the 
vicinity  of  a  third  harmonic  overtone  of  thickness  shear  mode 
in  a  SC-cui  quartz  strip. 

The  straight  crested  waves  are  spatially  one  dimensional,  but 
since  the  study  is  part  of  a  wider  study  of  third  overtone  SC- 
cut  quartz  resonators,  two-dimensional  finite  element  equations 
for  piezoelectric  plate  vibrations  are  derived  and  presented,  A 
one-dimensional  finite  element  program  is  developed  which 
could  be  modified  subsequently  to  handle  two-dimensional 
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problems.  The  mechanical  stiffness  of  the  finite  element  matrix 
equations  for  the  SC-cut  quartz  is  a  few  orders  of  magnitude 
larger  than  the  piezoelectric  stiffness,  hence,  the  piezoelectric 
coupling  is  weak,  and  the  piezoelectric  effects  could  be  mod¬ 
eled  using  a  perturbation  method.  A  new  perturbation  method 
for  piezoelectric  vibrations  is  derived  and  prop>oscd.  Results 
for  the  perturbation  method,  free  vibrations,  and  steady-state 
vibrations  of  a  partially  electroded  SC-cut  strip  are  shown. 

n.  Plate  Equations  for  High  Frequency 
Piezoelectric  Vibrations 

The  piezoelectric  plate  equations  for  high  frequency  vibra¬ 
tions  employed  in  the  current  finite  element  formulation  were 
derived  by  Lee,  Syngellakis,  and  Hou  [1],  The  components  of 
mechanical  displacement  and  electric  potential  are  expanded 
in  an  infinite  series  with  their  thickness-dependence  expressed 
by  trigonometric  functions: 

oo 

Ui(xi,X2,X3,f)  =  ^  uj  ^(xi,X3,i) 
n=0 

-  COS  ^(1  ““  0 

oo 

4>{.Xi,X2,Xi,i)  =  1,13,0 

n=0 

•cos  ^(1-0  (lb) 

where  ^  =  X2/6  is  a  nondimensional  thickness  coordinate 
of  the  plate.  The  terms  and  are  functions  of  spatial 
coordinates  xi  and  X3  only,  and  time  L  We  present  without 
formal  derivations  the  following  field  equations  from  [1],  The 
interested  reader  may  obtain  details  of  their  derivations  from 
the  same  reference.  Some  equations  are  given  in  a  matrix  form 
to  facilitate  the  development  of  finite  element  matrix  equations. 

A.  Stress  Equations  of  Motion 

J,(n)  _  ^  1  r(n)  ^  (1  +  (2a) 

where  the  terms  and  Tj,  are  n-th  order  stresses  and  are 
defined  as  the  following: 

i;W=  (2b) 

=  (2=) 
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A  high  frequency,  piezoelectric,  laminated  plate  theory  is  developed  and  presented  for  the 
purpose  of  modeling  and  analyzing  piezoelectric  thin-film  resonators  and  filters.  The  laminated 
plate  equations  are  extensions  of  anisotropic  composite  plate  theories  to  include  piezoelectric 
effects  and  capabilities  for  modeling  harmonic  overtones  of  thickness-shear  vibrations. 
Two-dimensional  equations  of  motion  for  piezoelectric  laminates  were  deduced  from  the 
three-dimensional  equations  of  linear  piezoelectricity  by  expanding  the  mechanical 
displacements  and  electric  potential  in  a  series  of  trigonometric  function,  and  obtaining  stress 
resultants  by  integrating  through  the  plate  thickness.  Relations  for  handling  the  mechanical  and 
electrical  effects  of  platings  on  the  top  and  bottom  surfaces  of  the  laminate  are  derived.  A  new 
matrix  method  of  correcting  the  cutoff  frequencies  is  presented.  This  matrix  method  could  also 
be  used  to  efficiently  correct  the  cutoff  frequencies  of  any  nth  order  plate  laminate  theories  which 
employ  Mindlin’s  form  of  polynomial  expansion  of  mechanical  displacements  and  electric 
potential  through  the  plate  thickness.  The  first  order  laminated  plate  theory  with  correction 
factors  for  cutoff  frequencies  and  slope  of  the  flexural  branch  at  large  frequencies  was  applied  to 
a  2-layer,  zinc  oxide-silicon  strip,  and  a  3-layer,  zinc  oxide-zinc  oxide-silicon  strip.  Open  circuit, 
dispersion  relations  were  generated  for  a  range  of  volume  fractions  and  compared  to  the  exact 
dispersion  relations.  Both  the  2-layer  and  3-layer  strip  show  similar  qualitative  comparison:  The 
present  theory  compares  fairly  well  with  the  exact  dispersion  relation  for  real  wave  numbers  and 
nonpropagating  (imaginary)  wave  numbers  which  are  smaller  than  0.5/.  The  extensional 
branch,  and  thickness-shear  branches  begin  to  deviate  from  the  exact  solution  when  the 
nondimensionalized  frequency  is  greater  than  one.  Consequently,  to  maintain  accuracy  of 
solutions  when  using  the  present  first  order  laminated  plate  theory,  one  should  limit  the 
calculation  of  resonant,  nondimensionalized  frequencies  to  less  than  1.1.  Results  for  the 
frequency  spectrum  in  the  vicinity  of  the  open  circuit,  fundamental  thickness-shear  frequency, 
and  modes  shapes  were  presented  for  the  2-layer  strip  with  fixed  edges,  and  a  volume  fraction  of 
silicon  equal  to  0.2.  The  technically  important  fundamental  thickness-shear  mode  is  found  to 
have  the  shear  component  strongly  coupled  with  extensional  component. 


I.  INTRODUCTION 

We  discuss  in  this  paper  the  development  and  applica¬ 
tion  of  a  high  frequency  piezoelectric,  laminated  plate  the¬ 
ory  intended  for  analyzing  and  modeling  piezoelectric  thin- 
film  resonators.  An  overall  view  of  the  thin-film  resonator 
technology  is  given  in  Refs.  1  and  2.  The  resonator  itself  is 
also  referred  to  by  other  names  such  as  bulk  wave  com¬ 
posite  resonator^  and  film  bulk  acoustic  resonator.^  Theo¬ 
retical  analysis  and  modeling  of  the  thin-film  resonator  was 
restricted  to  simple  one-dimensional  models  until 
Milsom’s^  two-dimensional  analysis  of  a  thickness- 
extensional  mode  resonator  employing  zinc  oxide  thin 
films  on  silicon;  and  later,  Tiersten  and  Stevens’s"^  analysis 
of  a  thickness-extensional  trapped  energy  resonator  with 
rectangular  electrodes  in  a  similar  device.  Piezoelectric 
plate  theories  for  low  frequency  flexural  and  torsional  vi¬ 
brations  of  piezoelectric,  laminated  polymer  plates  were 
studied  by  Ricketts^  and  Lee  and  Moon.^ 

The  thin-film  resonators  of  interest  are  excited  at  high 
frequencies  in  the  vicinity  of  the  fundamental  thickness 


modes  and  their  higher  harmonic  overtones.  The  analysis 
of  such  modes  of  vibrations  would  require  either  a  three- 
dimensional  theory,  or  a  high  frequency,  piezoelectric, 
laminated  plate  theory.  One  could  argue  that  a  piezoelec¬ 
tric  composite  could  always  be  modeled  using  three- 
dimensional  finite  elements  such  as  those  presented  by  Al- 
lik  and  Hughes’  and  Lerch.^  There  exist  commercial  finite 
element  software  packages  which  have  a  three-dimensional 
piezoelectric  element.  Such  a  route  for  high  frequency 
problems  is  extremely  intensive  in  terms  of  memory  re¬ 
quirements  and  numerical  computations.  The  difficulty  in 
identifying  the  resonant  frequency  of  interest  and  its  mode 
shape  increases  rapidly  with  the  frequency  of  vibration. 
The  problem  becomes  quite  intractable  for  higher  har¬ 
monic  overtones  of  thickness-shear  vibrations.  A  high  fre¬ 
quency  plate  theory  would,  besides  being  less  intensive  in 
terms  of  memory  requirements  and  numerical  computa¬ 
tions,  alleviate  part  of  this  problem  because  the  displace¬ 
ments  through  the  plate  thickness  are  assumed  known  and 
are  expressed  in  a  fixed  number  of  displacement  compo¬ 
nents  Hence,  other  than  its  theoretical  complexities 
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ABSTRACT 

Two  algorithms  for  eigenvalue  problems  in 
piezoelectric  finite  element  analyses  are  introduced.  The 
first  algorithm  involves  the  use  a  Lanezos  eigensolver, 
while  the  second  algorithm  uses  a  Rayleigh  quotient 
iteration  scheme.  In  both  solution  methods,  schemes  arc 
implemented  to  reduce  storage  requirements  and  solution 
time.  Also,  both  solution  methods  seek  to  preserve  the 
sparsity  structure  of  the  stiffness  matrix  to  realize  major 
savings  in  memory.  In  the  Lanezos  solution  method,  the 
structurjil  pattern  of  the  consistent  mass  matrix  is 
(ixj)loited  to  gain  savings  in  both  memory  and  solution 
time.  In  the  Rayleigh  quotient  iteration  method,  an 
algorithm  for  generating  good  initial  eigenpairs  is 
employed  to  improve  significantly  its  overall  convergence 
rate,  and  convergence  stability  in  the  regions  of  closely 
.sj)aced  eigenvalues  and  repeated  eigenvalues.  The  initial 
eigenvectors  are  obtained  by  interpolation  from  a  coarse 
me.sh.  In  order  for  this  iterative  method  to  be  effective,  an 
eigenvector  of  intercist  in  the  fine  mesh  muHt  rcHcmblc  an 
eigenvector  in  the  coarse  mesh.  Hence,  the  method  is 
effective  for  finding  the  set  of  eigenpairs  in  the  low 
frequency  range,  and  not  in  the  high  frequency  range 
where  the  eigenvectors  of  the  coarse  mesh  docs  not 
re.s(;ml)le  w(dl  tlieir  counterpurts  in  the  fine  mesh.  HeKultH 
of  example  problems  are  presented  to  show  the  savings  in 
solution  time  and  storage  requirements  of  the  proposed 
algorithms  when  compared  with  the  conventional 
algorithm  which  uses  static  condensation  to  remove  the 
electric  potential  degrees  of  freedom  from  the  piezoelectric 
eigenvalue  problem.  The  disadvantage  of  this  conventional 
scheme  is  that  the  static  condensation  destroys  the 
sparsity  structure  of  the  stiffness  matrix,  which  leads  to 
major  increases  in  memory  and  solution  time. 

INTRODUCTION 

The  analyBis  and  design  of  piezoelectric  resonator 
devices  present  a  challenging  set  of  attrib\itc8  to  the 
engineer:  (1)  The  geometry  and  boundary  conditions  of 
the  devices  are  complex,  (2)  the  resonating  element  is  a 
compoHilc,  conRisting  of  metal  electrodes  on  a  piezoelectric 
pliile  or  Hulrntrate,  and  in  thin  film  <leviccH,  the 


1051-0117/94/0000-1057  $4.00  ©  1994  IEEE 


piezoelectric  layer  is  grown  on  a  dielectric  layer,  and  (3) 
the  devices  usually  operate  at  higher  inodes  vibrations. 
Hence,  the  analysis  of  these  devices  demands  good 
analytical  models,  and  detailed  finite  element  solution  of 
such  analytical  models.  This  is  especially  relevant  when 
accurate  prediction  of  the  unwanted  modes  is  required. 

The  numerical  models  would  generate  a  set  of 
coupled  piezoelectric  eigenvalue  equations  of  the  following 
form: 


Afu  -f  “f  =  F, 

(1) 

(2) 

where  the  matrices  Af,  ATuu,  Ku(f>  and  arc  the  ma.ss, 

mechanical  stiffness,  piezoelectric  coupling,  and  dielectric 
matrices  respectively.  The  vectors  ta,  <f>,  F  and  Q  are 
respectively  the  mechanical  displacement,  electric 
potential,  mechanical  force  and  electric  charge  vectors.  A 
static  condcnsationfl]  of  the  potential  <f>  in  cqs.  (1)  and  (2) 
renders  a  single  equation  written  as 

Afu  +  ii:*u  =  F\ 

(3) 

where 

K* 

H) 

and 

P*  =  F  -  K.^K~lQ. 

(5) 

K*  is  the  condensed  clcctroclastic 

:  stiffness  matrix  and  F* 

is  the  corresponding  electromechanical  forcing  function. 
The  components  of  the  electric  potential  vectors  are 
recovered  from  cq.  (2)  by 

4>  =  k;;{q-k^.u).  (6) 

A  generalized  eigenvalue  problem  can  be  obtained  from  cq. 
(3)  by  assuming  harmonic  Holutions,  and  setting  F*  equal 
to  zero  to  yield 

(K*-u^M)u=0,  (7) 

where  u  is  the  nat\iral  frequency. 
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Abstract  Of  The  Dissertation 
Exact  And  Approximate  Analysis  Of  The  Propagation 
Of  Acoustic  Waves  In  Layered  Piezoelectric  Plates 
by  James  T.  Stewart 
Dissertation  Director: 

Professor  Yook-Kong  Yong 

Exact  and  approximate  analysis  of  the  propagation  of  acoustic  waves  in  layered 
piezoelectric  plates  is  performed  with  specific  application  to  the  analysis  of  thin  film  reso¬ 
nators.  The  exact  analysis  is  performed  using  the  transfer  matrix  method.  By  definition,  a 
transfer  matrix  is  a  linear  transformation  which  maps  a  specified  vector  of  field  quantities 
from  one  point  to  another  in  a  material.  By  constructing  transfer  matrices  for  each  layer 
and  enforcing  inter-layer  continuity  of  pertinent  field  quantities,  a  simple  matrix  multipli¬ 
cation  rule  results  which  can  be  used  to  map  these  quantities  from  the  plate’s  lower  free 
surface  to  its  upper  free  surface.  Hie  result  is  a  problem  which  is  in  terms  of  the  field 
quantities  evaluated  on  the  plate’s  free  surfaces  only,  which  is  no  more  difficult  to  solve 
than  it  is  for  a  homogeneous  plate.  The  transfer  matrix  analysis  is  directly  applicable  to 
problems  of  thickness  excitation  and  it  will  be  shown  that,  with  some  modifications,  this 
method  can  be  used  to  study  problems  of  lateral  excitation.  The  approximate  analysis  is 
performed  using  a  high  frequency  laminate  theory.  The  laminated  plate  equations  are 
developed  using  two  different  series  expansions  for  the  thickness  dependency  of  displace¬ 
ments  and  electric  potential.  A  general  correction  method  is  developed  which  improves 
the  accuracy  of  the  approximate  solutions.  In  the  first  analysis,  a  power  series  expansion  is 
used,  and  later  a  trigonometric  series  expansion  is  applied.  Cutoff  behavior  as  well  as  gen¬ 
eral  dispersion  characteristics  for  2anc  Oxide  on  Silicon  thin  film  resonators  are  studied 
using  the  approximate  analysis  and  compared  to  the  exact  solutions.  In  the  final  analysis,  a 
layered  piezoelectric  strip  of  finite  length  is  studied  under  conditions  of  free  vibration  with 
mechanically  clamped  ends. 
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ABSTRACT  OF  THE  DISSERTATION 


Numerical  Analysis  of  The  Third  Overtone  of  Thickness 
Shear  Vibrations  in  SC-cut  Quartz  Resonators 

by  Zlicn  Zhang  ,  Pli.D. 

Dissertation  Director:  Professor  Yook-kong  Yong 

Exact  and  approximate  dispersion  relations  up  to  the  third  overtone  of 
tliickness  shear  are  developed  for  an  infinite  SC- cut  quartz  piezoelectric  plate  with  the 
triclinic  symmetry  (all  36  quantities  are  present  in  the  elastic  tensor)  using  the 
three-dimensional  piezoelectricity  and  a  two-dimensional  piezoelectric,  high  frequency 
plate  theory,  respectively*  By  comparing  the  results,  the  approximate  theory  is  found 
to  provide  an  excellent  approximation  of  the  vibrational  behavior  of  SC-cut  quartz 
plate  resonators. 

A  realistic  and  detailed  finite  element  model  of  piezoelectric  plate  resonators, 
which  takes  into  consideration  of  the  effects  of  the  electrodes,  is  developed.  The  finite 
element  formulation  is  derived  based  on  the  general  N-th  order  plate  equations  using 
the  variational  principle.  The  finite  element  matrix  equations  from  the  tliird  order 
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plate  equations  are  used  in  the  analysis  of  the  third  overtone  of  thickness  shear.  The 
finite  element  model  developed  in  this  work  may  be  employed  in  high  frequency 
SC-cut  quartz  plate  resonator  analysis  to  obtain  characteristics  of  vibration  such  as 
resonant  frequencies  and  mode  shapes  wliich  would  otherwise  be  available 
experimentally  through  elaborate  and  costly  laboratory  work. 

Numerical  strategies  for  solving  very  large  scale  finite  element  model  of  SC-cut 
quartz  resonators  are  explored.  A  perturbation  technique  is  developed  to  account  for 
the  piezoelectric  effect  for  weakly  coupled  material.  This  method  saves  substantial 
memory  and  computational  time,  and  is  found  to  be  rather  accurate  for  SC-cut 
quartz.  In  the  calculations  of  the  purely  mechardcal  eigenpairs,  a  new  storage  scheme 
for  the  mass  matrix  is  proposed.  The  scheme  makes  use  of  the  special  structure  of  the 
mass  matrix  and  reduces  the  memory  requirement  by  more  than  90%  over  that  of  the 
envelope  storage  scheme  in  the  current  analysis.  The  substructuring  technique  is 
introduced  into  the  eigenvalue  calculations  to  further  diminish  the  memory  demands 
as  well  as  the  computational  time  without  any  loss  of  accuracy. 

Numerical  analysis  of  the  third  overtone  of  thickness  shear  of  SC-cut  quartz 
resonators  is  performed  using  the  finite  element  method.  Frequency  spectra,  resonant 
frequencies  and  mode  shapes  are  calculated  and  presented.  Vibrational  characteristics 
of  the  predominant  mode  is  studied  with  the  effects  of  the  electrodes. 


